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The energy density wake produced by a heavy quark moving through a strongly coupled N = 4 
supersymmetric Yang-Mills plasma is computed using gauge/string duality. 



Introduction. — The discovery that the quark-gluon 
plasma produced in heavy ion collisions at RHIC be- 
haves as a nearly ideal fluid [U [2] has prompted much 
interest in understanding the dynamics of strongly cou- 
pled non-Abelian plasmas. Gauge/string duality [3] H] 
allows one to compute many observables probing non- 
equilibrium dynamics of thermal J\f = 4 supersymmetric 
Yang-Mills (SYM) theory, including the rate of energy 
loss of a heavy quark moving through an SYM plasma 
[5]. (See also Refs. [7J HI EE EH] and references therein.) 
Energy transferred to the plasma from the moving quark 
will cause the energy density of the plasma, in the vicinity 
of the quark, to deviate from its equilibrium value. That 
is, the moving quark will create an energy density "wake" 
which moves with it through the plasma. The structure 
of this wake is of interest for studies of jet quenching and 
jet correlations in heavy ion collisions [HJ[I2]. We eval- 
uate this energy density perturbation, (AT 00 (x)), and 
display the resulting energy density wake in the case of 
subsonic, transsonic, and supersonic motion [20] . 

Gravitational description. — According to gauge/string 
duality, the addition of a massive quark to the J\f = 4 
SYM plasma is accomplished by embedding a D7 brane 
in the AdS-Schwarzschild (AdS-BH) geometry and then 
adding a string running from the D7 brane down to the 
black hole horizon [S] [T3]. The presence of the string 
perturbs the geometry via Einstein's equations. The be- 
havior of the metric perturbation near the AdS boundary 
encodes the change in the SYM stress-energy tensor. In 
the N c — > oo limit, the 5d gravitational constant becomes 
parametrically small and consequently the presence of 
the string acts as a small perturbation on the AdS-BH 
geometry. To obtain leading order results in Nf/N c , we 
write the full metric as g M „ = gffl + h^, where gffl is the 
metric of the AdS-BH geometry, and then linearize the 
Einstein equations in the perturbation h^. This gives 
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where A"£ is a second order linear differential operator, 
k 2 = 47T N 2 with L the AdS curvature radius, and t^ is 
the five dimensional stress-energy tensor of the trailing 
string solution [5j. The boundary value of the metric 
perturbation acts as a source for the SYM stress-energy 



tensor via the relation [2] 
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where Sb is the on-shell gravitational boundary action 
with h^ u the boundary value of the metric perturbation. 

We choose a coordinate system such that the metric of 
the AdS-BH geometry is 
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where f(u) = 1 — (u/uh) A - The event horizon is located 
at u = Uh, with T = (iruh)^ 1 the temperature of the 
SYM plasma. The boundary of the AdS-BH spacetime 
is at u=0. The geometry is translationally invariant in 
the four Minkowski space directions (t,x), so one may 
perform a spacetime Fourier transform and work with 
mode amplitudes h^ v (u\ uj, q). 

Gauge invariants. — The Ad J\f = A SYM stress-energy 
tensor is tracclcss and conserved, and consequently con- 
tains five independent degrees of freedom. This contrasts 
with the fifteen degrees of freedom contained in the met- 
ric perturbation h^ v . However not all of these degrees of 
freedom are physical. The linearized field equations are 
invariant under coordinate transformations — > af^+f*, 
where £ M is an arbitrary infinitesimal vector field. Under 
such transformations, the metric perturbation transforms 
as 



(4) 



where is the covariant derivative with respect to the 

background metric gjxJ. Consequently, five components 
of may be eliminated via gauge fixing. As in elec- 
tromagnetism, this does not completely fix the gauge. If 
one chooses to set h§ u = 0, where x$ = u denotes the 
AdS radial coordinate, then the residual gauge freedom 
allows one to eliminate five more components of h^ u on 
any single u = const, hypersurface. This limits the num- 
ber of independent physical degrees of freedom carried 
by h^v to five, matching that of the SYM stress tensor. 

Useful gauge invariants may be constructed out of 
linear combinations of the Fourier mode amplitudes 
h^(u;uj,q) [15], and classified according to their behav- 
ior under spatial rotations. There is one helicity zero 
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gauge invariant linear combination, and a pair each of 
helicity one and two invariants [21] . Rotation plus gauge 
invariance implies that these five invariants satisfy de- 
coupled equations of motion. For later convenience, let 
H^ v = {u 2 / L 2 ) h^v A short exercise using Eq. Q shows 
that 

Z(u;uJ,q) = q 2 H o + 2ujq l H t0 + (uj 2 /q 2 )q l q'H tJ 

+ \ [2-f(u)-u 2 /q 2 ] tfHu-jqiHij) (5) 

is the helicity zero gauge invariant linear combination of 
metric perturbations, unique up to multiplication by an 
arbitrary function of u. (Repeated Minkowski spatial in- 
dices i,j = 1,2,3 are implicitly summed.) It is straight- 
forward (but tedious) to work out the equation of motion 
for Z from the linearized field equations. Doing so, and 
inserting the explicit form of the trailing-string gravita- 
tional stress tensor (also computed in Ref. [8]) yields the 
second order ODE 



Z" + A{u) Z' + B{u) Z = S(u) , 



where 



A{u) 

B{u) 
S{u) 



uf 
f 



24 (q 2 f- 



q 2 (u/'-6/) + 6w 2 



10 

T 



32 q 2 u 6 ul 



_ n 5 



K 2 V\ q 2 (v 2 - 
6ttL 3 



q 2 (uf'-6f) 
2) -3uj 2 



6uj 2 



(6) 

(7) 
(8) 



q 2 Vl 

.2 



2„5 



UJUtU 



9(g 



2 -u 2 ) 2 ul 



f(fq 2 + 2q 2 - 
x 2tt(5(w-v • q) e- iuJ x ' t " 



(9) 



Here A is the t' Hooft coupling, v is the quark velocity, 
and 
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is the trailing string profile [5] . 

In a gauge in which h^ a —0 for all a, the change in the 
energy density is given by [T3] 
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where is the quartic term in the expansion of 

H^ v = (u 2 /L 2 ) at the boundary [22]. A simple con- 
nection between Z and the energy density may be found 
by considering the behavior of Z and H^ v near u = 0. 
Substituting a power series expansion into the linearized 
equations of motion ([!]) , one finds that near the boundary 
Z and H^ v have the forms 



Z(u) = Z( 3 ) U 3 + Z( 4 ) u A + 



(12) 
(13) 



Moreover, the three combinations q % q 2 H^f\ and 

q l q^H < lf > are all determined by Hqq . Substituting the 
coefficients Hjn) into Eq. (I5l and solving for yields 
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where 



A 



_ 4V\ iu; [g 2 (5v 2 +l) - 3 (v 2 +l) oj 2 
~ 8ttL 3 q 2 u 2 Vl^v 2 
x (2-7r) (5(cj— v • q). 



(15) 



The temperature dependent perturbation in the energy 
density (which we focus on below) is therefore given by 
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where AZ^ is the difference between Z^ evaluated at 
temperature T and zero temperature. 

Asymptotics and numerics. —We solve Eq. ^ with 
a Green's function G(u, u') constructed out of homoge- 
neous solutions. 



G{u,u')=g < {u < )g > {u > )/W{u'), 



(17) 



where W(u) is the Wronskian of g< and g > . The appro- 
priate homogeneous solutions are dictated by the bound- 
ary conditions. The differential operator in (|6| has sin- 
gular points at u — and u — with exponents and 4, 
and iiwwh/4, respectively. Vanishing of the metric per- 
turbation near the boundary requires that g < (u) ~ u 4 as 
u — > 0, while the requirement that the black hole not ra- 
diate [16 implies that g>(u) ~ (u~Uh)^ luJUh ^ 4 near the 
horizon. The overall normalization of <?< may be fixed 
by requiring lim„^o g<(u) /u A = 1. Zero temperature so- 
lutions to Eq. (|6| are easily found and involve modified 
Bessel functions. A short exercise leads to 



AZ, 



(4) 



lh du {?Ms{u) 



< 

I W(u) 



T . . q 2 ~LU 2 T ^ ( a , 



^)*<->}, da) 



where J = (/ — uf')/ f 2 is a Jacobian factor and So(m) 
is the source 5*(u) evaluated at T = 0. 

For very large or small momneta (compared to T), one 
may find explicit asymptotic expressions for the homo- 
geneous solutions and derive the asymptotic behavior of 
£(o->, q). The large momentum limit is physically identi- 
cal to the low temperature limit, so expanding Eq. |6]) 
about Uh = oo allows one to extract the large q asymp- 
totics. For a fixed ratio r = ui/q, the linear differential 
operator on the left side of Eq. (|6| differs from its T=0 
limit only by terms of relative order T^/q 4 , while the 
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FIG. 1: Plot of |x| £(x)/(T_VT) for v = 1/4, with the zero FIG. 2: Plot of jxj £(x)/(T 3 VA) for v = 3/4, with the T=0 
temperature and near zone (20| contributions removed. Note an d near zone (poj) contributions removed. A Mach cone is 
the absence of structure in the region |x| > 1/wT. clearly visible, with an opening half-angle 6 w 50°. 



source S equals So (l + ioju 3 /3u 2 ) up to 0(T 4 /q 4 ) cor- 
rections. Denning for convenience 



s(w,q) 
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we find the large momentum asymptotic behavior |23j 

eM = sM-^ — If- \ — '—L, (20) 

up to relative corrections suppressed by 1/q 2 . For small 
momentum, we expand Eq. (|6j) in powers of q with the 
ratio uj/q fixed. The linear differential operator on the 
left side of Eq. ^ has a smooth q — ► limit with 0(q 2 ) 
corrections, but relative O(q) terms appear in the source 
S and the Green's function (via the incoming boundary 
condition at the horizon). As q — > 0, we find 
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up to relative corrections suppressed by q 2 . 

To compute the energy density £ (x) for a given quark 
speed i!=|v|, we use the cylindrical symmetry of the 
source ^ to reduce the spacetime Fourier transform 
to a two dimensional integral over qu = v • q and 
q± = |v x q|, which is performed numerically. For each 
value of (/|| and q±, the Fourier amplitude £(vq\\, qti, qx), 
as given by Eq. (16 1, is evaluated by numerically inte- 
grating the homogeneous differential equation ^ (with- 
out source) outward from the horizon to find g > (u), and 



then evaluating numerically the radial integral (18 1 to 



find AZu\(vq\\,q\\,q±.). (The u dependence of the Wron- 
skian can be computed analytically.) 

Results and discussion. — For small distances d = 
|x— vt| <C 1/T away from the moving quark, the dom- 
inant contributions to the energy density come from mo- 
menta g > T. The leading short distance behavior is 
temperature independent; T=0 conformal invariance im- 
plies that the energy density scales like 1/d 4 . The first 
temperature dependent near zone contribution to the en- 
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FIG. 3: Plot of |x|£(x)/(T 3 VA) for v = l/VS, with the 
T=0 and near zone ( 20 1 contributions removed. A planar 



shock front perpendicular to the quark velocity is evident. 



sition space energy density which scales like T 2 / d 2 in the 
vicinity of the quark. 

Figures [l||3] show real space plots of £ (x) for quark ve- 
locities v = 1/4, 3/4, and 1/V3, respectively, with the 



ergy density comes from the term ( 20 1 . This yields a po- 



near zone contribution ( 20 1 removed (in order to high- 

M 



light the intermediate and far zone structure) |24j . In 
these plots the quark mass is infinite and the quark's lo- 
cation, at the time shown, is x — 0. Since Af = 4 SYM is 
a conformal theory, the speed of sound is 1/ %/3- Hence 
Fig. 1 shows subsonic motion, Fig. 2 shows supersonic 
motion, and Fig. 3 is precisely at the speed of sound. 
For all three velocities we observe a net surplus of en- 
ergy in front of the quark and a net deficit behind the 
quark. This may naturally be interpreted as plasma be- 
ing pushed and displaced by the quark, just like the water 
displacement produced by a moving boat. 

The most striking feature in these plots is the forma- 
tion of a conical energy wake, or sonic boom, for velocities 
v > v s = l/\/3. A textbook constructive interference ar- 
gument shows that a projectile moving supersonically, in 
any fluid, should produce a Mach cone with an opening 
half-angle given by sin# = v s /v (where tan# = — xj_/x||). 
For v — 3/4, this is 50.3°. For the transsonic v—v s case 
shown in Fig. [3] we see an energy wake along the plane 
front £|| = 0, while for v = 3/4 the wake is concentrated, 
as expected, along a 50° cone. From Fig. [2] one may see 
that the shock wave has a width w 10/7rT and broadens 
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with increasing distance. This behavior is to be expected 
in a viscous fluid where sound waves are damped and 
diffusive. Furthermore, the amplitude of the shock front 
decreases slightly faster than 1 jd. If the shock front did 
not broaden with increasing distance, then conservation 
of energy would imply a 1/d decrease of the intensity of 
the shock wave in the far zone. 

It is instructive to compare the long wavelength limit 



(21 1 of the energy density with the behavior predicted by 
linearized hydrodynamics, in which the energy density 
perturbation £ satisfies the diffusive wave equation 



7 V 2 <9 t 



£ = p. 



(22) 



Here p is an effective source which depends on dynamics 
in the near zone, and 7 = is the sound attenuation 

constant. (See also Ref. [17].) In strongly coupled N = 4 
SYM, 7 = (37rT) -1 [IS]. Fourier transforming and ex- 
panding in powers of momentum (at fixed r = uj/q) with 



p= qpi+q p2 



£ = - 



3pi 
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, one has 

(l-3r 2 ) 2 
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Comparing the above with the sound poles in Eq. (211, we 



see that our result for the far zone energy density agrees 
with linear hydrodynamics provided 7 has the expected 
value of (37rT) _1 and we identify 



p(x,t) = - 



(1+« 2 )9*<5 3 (x-v<), 



(24) 



up to higher derivative corrections. This may be com- 
pared with the source obtained from the energy momen- 
tum conservation equation d^T^ = F u , where — 
<5 3 (x— vt)y/l— v 2 and / M is the external force (or mi- 
nus the drag force) acting on the quark. Using the equa- 
tions of linearized hydrodynamics, one may calculate p in 
terms of F v . To leading order (in derivatives) one finds 



V -F 



d t F°. 
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Inserting the drag force computed in Ref. 5 precisely re- 
produces the result ( 24 1 . Unsurprisingly, higher deriva- 



tive terms in the source are not reproduced by linear 
hydrodynamics. 

Natural extensions of this work include the calculation 
of energy flux and the inclusion of a non-zero chemical po- 
tential. As this work was nearing completion, we learned 
of a similar study [19] of the energy density of a moving 
quark. We thank the authors of Ref. [19], as well as A. 
Karch, CP. Herzog, and D.T. Son, for useful discussions. 
This work was supported in part by the U.S. Department 
of Energy under Grant No. DE-FG02-96ER40956. 
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